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ABSTRACT 

We have used the Sydney University Stellar Interferometer (SUSI) to measure the angular 
diameter of the F9 V star (3 Virginis. After correcting for limb darkening and combining with 
the revised Hipparcos parallax, we derive a radius of 1.703 ± 0.022 i?0 (1.3 per cent). We 
have also calculated the bolometric flux from published measurements which, combined with 
the angular diameter, implies an effective temperature of 6059 ± 49 K (0.8 per cent). We also 
derived the luminosity of /3 Vir to be L = 3.51 ± 0.08 ^ (2.1 per cent). Solar-like oscil- 
lations were measured in this star by Carri er et al.l (|200?) and using their value for the large 
frequency separation yields the mean stellar density with an uncertainty of about 2 per cent. 
Our constraints on the fundamental parameters of Vir will be important to test theoretical 
models of this star and its oscillations. 

Key words: stars: individual: (3 Vir - stars: fundamental parameters - techniques: interfero- 
metric 



1 INTRODUCTION 

Accurate measurements of fundamental stellar parameters - par- 
ticularly radius, luminosity and mass - are essential for critically 
testing the current generation of stellar models. Some of these pa- 
rameters can be constrained using well-established techniques such 
as spectrophotometry, parallax and optical interferometry. In the 
last decade, observing stellar oscillation frequencies (asteroseis- 
mology) has been shown to be a powerful tool to infer the inter- 
nal structure of stars (e.g. iBrown & Gillilandlll994) . In particular, 
the mean stellar density can be estimated from the large frequency 
separation between consecutive radial overtones. 

The links between asteroseismology and interferometry have 
been comprehensively reviewed by Cunhaetal. (2007). Recently, 
the complementary roles played by these two techni ques in the in- 
vestig ation of solar-type stars have been described bv lCreevev et al.l 
( I2007I) . who showed the mass of a single star can be determined to 
a precision (l-a) of better than 2 per cent by combining a good 
radius measurement with the density inferred from th e large sep- 
aration . Indeed, this precision was almost obtained bv lNorth et al.l 
when they constrained the mass of (3 Hyi to a precision of 
2.8 per cent using re sults from optical in terferometry and astero- 
seismo logy (see also lKieldsen et aL| l2008). Moreover. iNorth et al.l 
( l2007t) further constrained the fundamental parameters of /3 Hyi 
with values for the stellar radius, luminosity, effective temperature 
and surface gravity. 

The star /3 Virginis (HR 4540, HD 102870, HIP 57757) has 
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spectral type F9 V jMorgan & Keenarilll973h and is slightly metal- 
rich, with [Fe/H] = 0.14. It is bright {V = 3.60) and has a w ell- 
determined parallax. Models bv lSggenberger & Carried ( 1200 6) us- 
ing the Geneva evolution code, including rotation and atomic dif- 
fusion, indicate that the star has a mass of 1.2-1.3 Mq and is close 
to, or just past, the end of main-sequence core-hydrogen burn- 
ing. Evidenc e for s olar-like oscillations in /? Vir was reported by 
iMartic et al.1 ( i2004l) . in the for m of excess power at frequencies 
around 1.7 mHz. Subsequentlv. ICarrier et al. |(|2005|) confirmed the 
presence of oscillations and identified a series of regularly spaced 
modes with a large frequency separation of Av — 72.1 /iHz. 

In this paper we present the first interferometric measurement 
of the angular diameter of (3 Vir (Section[2}- We also calculate the 
bolometric flux of the star using data from the literature (Section[3](. 
These, together with the parallax from Hipparcos and the density 
from asteroseismology, allow us to constrain the fundamental stel- 
lar parameters of f3 Vir (effective temperature, radius, luminosity 
and mass; Section|4ll. 



2 ANGULAR DIAMETER 

2.1 Interferometric Observations with SUSI 

Measurements of the squared visibility (the normalised 
squared-modulus of the complex visibility) were made on a to- 
tal of 1 2 nights using the Sydney University Stellar Interferometer 
(SUSI: lOavis et aUll999h . We recorded interference fringes with 
the red-table beam-combination system using a filter with (nom- 
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Table 1. Parameters of calibrator stars. 



HR 


Name 


Spectral 


V 


UD Diameter 


Separation 






Type 




(mas) 


from /3 Vir 


4368 


(f> Leo 


A71V 


4.47 


0.58 ± 0.08 


10?08 


4386 


cr Leo 


B9V 


4.04 


0.42 ±0.06 


8? 52 


4515 


5 Vir 


A4V 


4.84 


0.45 ±0.06 


6° 63 



inal) centre wavelength 700 nm and full-width at half-maximum 
80 nm. Th is beam-combination system has been described in 
lOavis et a l. (2007a) along with details of the standard SUSI ob- 
serving, data reduction and calibration procedures. 

The parameters of the calibrator stars are given in Table [T] 
There are no measured angular diameters for these stars, and the 
values given here were estimated from {B — V)o colours by in- 
terpolating measurements of similar stars with the N arrabri Stellar 

1 19741) and the Mark 



Intensity Interferometer (Hanburv B rown et al 



III Optical Interferometer (iMozurkewich et al 



20031) . Corrections 



for li mb-darkening were done using the results of iDavis et al] 

jiooo ). 

The journal of observations is given in Table |2] We obtained 
a total of 77 estimations of V^. It should be noted that, for consis- 
tency, the uncertainties in the angular diameters of the calibrating 
stars were included in the V'^ uncertainty calculation, even though 
they have a negligibly small effect when compared to the measure- 
ment uncertainty. 



2.2 Analysis of Visibilities 

In the simplest case, the brightness distribution of a star can be 
modeled as a disc of uniform irradiance with angular diameter Oijo . 
The theoretical response of a two-aperture interferometer to such a 
model is given by 



\V\' 



2Ji(7rB6»uD/A) 



TTBevu/X 



(1) 



where B is the projected baseline, A is the observing wavelength 
and Ji is a first order Bessel function. However, real stars are limb- 
darkened and corrections are needed to fin d the 'true' angular di- 
amete r. These corrections can be found in lOavis. Tango & Boothl 
1 I2OOOI) and are small for stars that have a compact atmosphere. 

To account for any systematic effects arising from the slightly 
differing spectral types of (3 Vir and the calibrators, the presence of 
a close companion, or other causes, an additional parameter. A, is 
included in the uniform disc model. Equation ([TJ becomes 



,2Ji(7rB6»uD/A) 



ttBOud/X 



(2) 



The interferometer response given in equation ^ is, 
strictly sp eaking, only valid for monochromatic observations. 
iTango & Davis ( 2002) have analysed the effect of a wide observing 
bandpass on interferometric angular diameters and found it to be 
insignificant provided the interferometer's coherent field-of-view 
is larger than the angular extent of the source. The coherent field- 
of-view of SUSI during observations was found to be greater than 
5.8 mas, hence greater than the extent of f3 Vir (see below). There- 
fore, bandwidth smearing can be considered negligible. However, 
SUSI's effective wavele ngth when observin g an F9 star is approxi- 
mately 696.4 ± 2.0 nm jPavis et al.ll2007ij) . The fit to equation ID 
was completed with this effective wavelength. 
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Figure 1. Nightly weighted-mean measures with the fitted uniform disc 
model overlaid. Inset: all measures at the (nominal) 80 m SUSI baseline. 



An implementation of the Levenberg-Marquardt method was 
used to fit equation {2) to all measures of to estimate 6'ud and 
A. The diagonal elements of the covariance matrix (which are cal- 
culated as part of the minimisation) were used to derive formal 
uncertainties in the parameter estimation. We used Markov chain 
Monte Carlo (MCMC) simulations to verify the formal uncertain- 
ties because equation l|2j is non-linear and the squared visibility 
measurement errors (derived by the SUSI reduction pipeline) may 
not strictly conform to a normal distribution. A MCMC simula- 
tion, implemented with a Metropolis-Hastings algorithm, involves 
a likelihood-based random walk in parameter space whereby the 
full marginal posteri or probability d ensity function is estimated for 
each parameter (see lGregorvll200j for an introduction to MCMC). 
An advantage of this method is that prior information (and the ef- 
fect on parameter uncertainties) can be included in the simulation. 
In the case of /3 Vir, the effect of the observing wavelength uncer- 
tainty was included in the simulation by the adoption of an appro- 
priate likelihood function (see below). 

The reduced of the original fit was 2.08, implying that the 
squared- visibility measurement uncertainties were underestimated. 
We have therefore scaled the measurement uncertainties by \/2.08 
to obtain a reduced of unity. The formal uncertainties were veri- 
fied with a series of MCMC simulations, each comprising 10^ iter- 
ations. For the effective wavelength we adopted a Gaussian likeli- 
hood function with mean 696 nm and standard deviation 2 nm. The 
probability density functions produced Gaussian distributions with 
standard deviation of similar values to the formal uncertainties. 
We therefore obtain (with the associated la formal uncertainty) 
SuD = 1.386 ± 0.017 mas and A = 0.985 ± 0.006. The data are 
shown in Fig. [T] with the fitted uniform-disc model overlaid. 

The fitted value of A differs slightly from unity, indicating 
there is either a close companion or some instrumental effect caus- 
ing a small reduction in the intercept. While [3 Vir does have 
known optical companions, these stars are too distant and too faint 
to affect the SUSI observations. Further more, /3 Vir appears in 
planet search programs (see for example Iwittenmver et alJl2006h 
and there is only a very limited parameter space where a further 
companion could reside without having already been detected. We 
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Table 2. Summary of observational data. The night of the observation is given in columns 1 and 2 as a calendar date and a mean 
MJD. The nominal and mean projected baseline in units of metres is given in columns 3 and 4, respectively. The weighted-mean 
squared visibility, associated error and number of observations during a night are given in the last three columns. 



Date 


MJD 


Nominal 
Baseline 


Mean Projected 
Baseline 


Calibrators 


y2 


cr 


#y2 


2007 February 14 


54145.70 


40 


34.05 


cr Leo, (; 


b Leo, g Vir 


0.684 


0.030 


3 


2007 February 15 


54146.66 


60 


51.28 


(J Leo, (; 


b Leo 


0.567 


0.023 


8 


2007 March 09 


54168.58 


5 


4.28 


fT Leo, (; 


b Leo 


1.006 


0.037 


3 


2007 March 10 


54169.58 


60 


51.35 


cr Leo, (, 


b Leo 


0.524 


0.029 


5 


2007 March 1 1 


54170.62 


15 


12.74 


(7 Leo, (; 


b Leo 


0.893 


0.018 


5 


2007 March 13 


54172.63 


5 


4.26 


cr Leo, (; 


b Leo 


0.956 


0.015 


7 


2007 March 14 


54173.56 


40 


34.47 


cr Leo, (; 


b Leo 


0.706 


0.015 


6 


2007 March 14 


54173.68 


15 


13.25 


cr Leo, (; 


b Leo 


0.951 


0.063 


2 


2007 March 15 


54174.59 


80 


68.16 


cr Leo, (; 


b Leo 


0.264 


0.010 


8 


2007 April 18 


54108.47 


5 


4.31 


cr Leo, (; 


b Leo 


0.998 


0.014 


10 


2007 April 20 


54110.49 


20 


17.04 


cr Leo, <; 


b Leo 


0.957 


0.032 


4 


2007 April 22 


54112.49 


80 


68.01 


cr Leo, (; 


b Leo 


0.339 


0.014 


6 


2007 May 27 


54147.39 


80 


68.32 


cr Leo, (; 


b Leo 


0.304 


0.009 


10 



also cannot explain the value of A in terms of a systematic reduc- 
tion of V'^ due to the slightly differing spectral types of target and 
calibrator stars. We believe it is more likely that the observations 
were affected by a small unknown systematic error or a binary 
companion of unusual characteristics, than that they were subject 
to unlucky large random errors. Therefore, we retain the original 
uniform disc angular diameter, Otju — 1.386 ± 0.017 mas, fitted 
with A as a free parameter. 

The correction for the effects of limb-darken ing was deter- 
mined to be Sld = 1.046SuD, using iDavis et alj feoOO ) and the 
following parameters: Tch = 6059 K, log g = 4.12 and [Fe/H] = 
0.14. With the exception of [Fe/H], these values are the result of an 
iterative procedure whereby initial values were refined by the sub- 
sequent determination of fundamental parameters in Section|4] The 
starting values for TcS and logg, and the fixed value for [Fe/H], 
were the mean values from those found in the literature, and are 
given in Table [3] The refinement of TcB and log g only caused a 
small change in the fifth significant figure of the limb-darkening 
correction, for which we estimate an uncertainty 0.002. Therefore 
the limb-darkened diameter is 6ld = 1.450 ± 0.018 mas and car- 
ries th e caveat that we are assuming that the Kurucz models, upon 
which iDavis et aL I i2000l) based their work, are accurate for this 
star. 



3 BOLOMETRIC FLUX 

We have determined the bolometric flux, /, based on flux- 
calibrated photometry and spectrophotometry from the litera- 
ture in combination w ith MARCS model stellar atmospheres 
(iGustafsson et al.l2003h . It has been determined by summing the in- 
tegrated fluxes over four spectral ranges, namely < 0.33 /im, 0.33- 
0.86 fim, 0.86-2.2 /im and 2.2-20 nm. The flux beyond 20 nm is 
< 0.01 per cent for a 6150 K black body and will be similar for 
the stellar flux distribution and can be ignored. It is noted that the 
(B — V) and (U — B) colours of (3 Vir are consistent with those of 
an u nreddened star o f its spectral classification. 

lKiehlinjh987i) has published spectrophotometry of (3 Vir for 
the wavelength range 325-865 nm. The observations were made at 
equal intervals of I nm with a resolution of 1 nm. The published 
spectral energy distribution is averaged over bandpasses 5 nm wide 



Table 3. Literature values for parameters of jS Vir. 



Source 


Toff (K) 


Method* 


logs 


[Fe/H] 


1 


6150 ± 100 


a 


4.2 ±0.1 


0.10 


2 


6190 ±80 


b 


4.20 


0.13 


3 


6176 


b 


4.14 


0.13 


4 


6127 ±55 


c 






5 


6124 ±31 


c 






6 


6124t 




4.24t 


o.igt 


7 


6068 ± 70 


b 


4.09 ±0.1 


0.13 ±0.1 


8 


6140 


d 


4.09 ± 0.08 


0.15 


9 


6055 ± 48 


e 






10 


6124 


f 




0.13 ±0.10 


11 


6076 ± 119 


b 


4.142 


0.13 


Mean 


6123 




4.16 


0.14 



So urce referen c es: (1) iThevenin. Vauclair & Vauclaiij il986l): 
Q) iBalachandranI h990l): (3) lEdvardsson et alj h993l): (4) Idi Benedetto! 
(1998); (5) Blackwell & Lynas-Gray' ("1998'): (6) Malagnini et al.' ('2OO0|); 
(7) Chenetal (2001); (8) Gray et al. (2001); (9) Kovtyukh et al. ( 20031) : 
(lO) IMorel & Micelal <2004l) : (1 D iAllende Prieto et al.l <2004l) . 
t; mean of listed values from elsewhere, many identical; 
#: method used for T^ff determination: (a) Fit to H7 profile; (b) From h 
and (6 — y) calibration; (c) IRFM method; (d) Fits of model spectra to 
spectra and uby; (e) From line depth ratios; (f) From iron line excitation 
and ionization equilibrium. 

and is tabulated e very 5 rmi, bas ed on the spectrophotometric cali- 
bration of Vega bv lHavesI fl98^ . This has been converted to a flux 
distribution using a value for the flux from Vega at 550.0 nm of 
3.56xI0~" Wm~^ nm~^ jMegessiei|[T99^ . In the case of 5 CMa 
dPavis et al."2007bY Kiehling's spectrophotometry was compared 
with that of Davis & Webb (1974) and the calibrated flux distribu- 
tions were found to be in excellent agreement, with an RMS dif- 
ference computed from the wavelengths in common of < 1.1 per 
cent and with no systematic differences over the wavelength range 
in common (330-808 van) . The fl ux distribution in the MILES li- 
brary of empirical spectra jSanc hez-Blazquez et al. 2006) for P Vir 
has also been considered and has been flux-calibrated in the same 
way as the Kiehling data. The wavelength coverage of the MILES 
flux distribution is 355-740 nm, less than the Kiehling range of 
325-865 nm. It is tabulated at 0.9 nm intervals with a resolution 
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of 0.23 nm, compared with the data by Kiehhng, which were aver- 
aged and tabulated over 5 nm intervals. The two distributions are in 
excellent agreement except for two apparently discrepant points in 
the Kiehling distribution at 760 nm and 765 nm. The two flux distri- 
butions were integrated for the common wavelength range of 355- 
740 nm, omitting the discrepant points, and the integrated fluxes 
agree to within 1 per cent. The Kiehling flux distribution covers a 
greater wavelength range and extends to the ultraviolet data at the 
short wavelength end and, for these reasons, it has been used to 
determine the integrated visual flux. 

The uncertainty in the integrated flux has been estimated by 
combining the uncertainty in the Megessier ( 1995) flux cali bration 
of 0.7 per cent, the uncertainty in the Vega calibration by iHavej 
l ll985h o f 1.5 per c e nt and the uncertainty in the relative flux distri- 
bution o f lKiehlind dlgSTl) . which is estimated to be ~1.5 per cent. 
This latter figure is based on the previous experience with 5 CMa 
and the good agreement between the Kiehling and MILES flux dis- 
tributions for /3 Vir. The resultant uncertainty is ±2.2 per cent and 
the integrated Kiehling flux for the wavelength range 0.33-0.86 /im 
is (5.73 ± 0.13) X 10-1° Wm-^ 

There ar e four flux values f or the ultraviolet from the 
TDl satellite jThompson et alj|l978i) at 156.5, 196.5, 236.5 and 
274.0 nm and these have been downloade d from SIMBAD. They 
have been plotted with the lKiehlind ( Il987l ) calibrated fluxes at 325 
and 330 nm and a smooth curve drawn through the six flux points. 
The flux shortward of 330 nm is only a small fraction of the total 
flux (<4 per cent), justifying this simple approach. The area un- 
der the curve shortward of 330 nm has been integrated to give the 
ultraviolet flux equal to (3.5 ± 0.2) x 10-^^ Wm"^ where the un- 
certainty is conservatively based on the published uncertainties in 
the ultraviolet fluxes. 

There are few infrared measurements for (3 Vir and it was 
necessary to interpolate and extrap olate them with the aid of a 
MARCS model stellar atmosphere jGustafsson etai] 1200 3!). The 
calibrated observational photometric fluxes are listed in Table [4] 
together with their sources and the references used for their cali- 
bration, and the fluxes are shown in Fig. (2] Initially the flux dis- 
tributions of MARCS model atmospheres for 6000 K and 6250 K, 
both for log g — 4.0, were scaled to fit the observational data in the 
wavelength range 0.6-2.2 /im. The observational data included the 
Kiehling fluxes from 0.6-0.865 /im and the fluxes for the RIJHK 
photometric bands listed in Table|4] The 6000 K model gave a good 
fit but the 6250 K model clearly did not. The scaled fluxes for the 
6000 K model were integrated in two ranges, 0.86-2.2 /im and 2.2- 
20 ^tm. The resulting integrated fluxes were added to the fluxes for 
<0.33 /xm and 0.33-0.86 /im determined above and, in combina- 
tion with the limb-darkened angular diameter, gave an effective 
temperature of -^6050 K. 

The fluxes for an effective temperature of 6050 K were inter- 
polated from the 6000 K and 6250 K models and the fitting pro- 
cedure repeated. The models have fluxes tabulated at intervals of 
0.005 per cent in wavelength, which results in very large plot files. 
For diagrammatic purposes the fitting was therefore also carried out 
with different flux averages — 0.025, 0.05 and 1.0 per cent intervals 
in wavelength. No significant difference was found in either the 
scaling factor for a fit or in the integrated fluxes for the different 
flux averages. The fit for the interpolated fluxes for 6050 K aver- 
aged over 0.025 per cent wavelength intervals is shown in Fig. [2] 
Uncertainties in the integrated fluxes have been derived by com- 
bining the estimated uncertainty in the fit to the observational data 
with the uncertainty in the observational data and in particular that 
in the Kiehling data involved in the fit. 



Table 4. Calibrated photometric IR fluxes for /3 Vir 



Band A^ff 
(pm) 



Flux 

(10-l2Wm-2 



Source Calibration 



/tm 



R 


0.641 


105.8 


1 


a 


I 


0.798 


71.6 


1 


a 


R 


0.70 


99.0 


2 


b 


I 


0.90 


61.0 


2 


b 


J 


1.25 


28.0 


2 


c 


K 


2.20 


4.8 


2 


c 


J 


1.2790 


27.05 


3 


d 


H 


1.6483 


12.54 


3 


d 


K 


2.1869 


5.01 


3 


d 


12/jm 


12 


0.0104 


4 


e 



So urce references: ( 1) ICousinj il980l): (2) Ijohnson et al] h966l) ; 
Ol lAlonso et al.Hl994 : r4lllRAS Team' ('1988'). 

Ca libration referenc es: (a)lBessell et al. I.1998h: (b) Ijohnsoril h966l) : 
(c^ lMegessieil Jl995h : fdl Alonso et al. 1 19941: fe^ llRAS Teanj jl988h re- 
duced by 4.1 per cent ICohen et al..l996.1 . 
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1.00 
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Figure 2. The flux distribution for /3 Vir for the wavelength range 0.5- 
2.5 /im, with data from Table|4] T he symbo l s have been plotted ov ersize for 
clarity. Key: small filled circles - lKiehlin3 (l987|) ; open squares - ICousinsl 
(1980) VRI photome try; filled squares - Johnson et al.l (196g) RI \iho- 
tometr y; open c ircles - Ijohnson et afl jl96(: ) JK photometry; filled circles 
- Alo nso et al ] il994 JHK photometry; grey line - fluxes for the fitted 
MARCS model atmosphere averaged over 0.025 per cent wavelength inter- 
vals. 



The integrated fluxes for the fitted 6050 K flux distribution are 
listed in Table [5] for the four wavelength bands considered. The 
uncertainties have been derived by combining the estimated uncer- 
tainty in the fit to the observational data with the uncertainty in 
the observational data and in particular that in the Kiehling data 
involved in the fit. 

Figure|3]shows an assembly of the observational flux data with 
the interpolated flux curve for a 6050 K model atmosphere fitted to 
the observational data as shown in Fig. |2] Figure [3] includes mid- 
infrared photometry from Spitzer by Trilling et al. (2008) and we 
note that the data appear to rule out their suggestion of a possi- 
ble excess from a debris disk. The 12 /tm IRAS Point Source flux 
does lie above the model curve but drawing a smooth curve from 
the flux for the K band through the IRAS point and integrating it 
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Figure 3. The flux distribution for /3 Vir from the ultraviolet to the mid- 
infrared. Symbols have the same meanings as Fig.|2] with the addition of: 
filled t riangles - IRAS P oint Source Fluxes; and open triangles - Spitzer 
fluxes iTrilling et ailEoOSl) . The dashed grey line is an extension of the 
MARCS model beyond 20 /im using a slope of A"'*. 



Table 5. Integrated fluxes for /3 Vir in each spectral band and the resulting 
bolometric flux (/). 



Wavelength 


Flux 


Range 


(lO-iOwm-2) 


(fim) 




0-0.33 


0.35 ±0.02 


0.33-0.86 


5.73 ±0.13 


0.86-2.2 


2.97 ±0.09 


2.2-20 


0.39 ±0.02 


>20 


negligible 


/ 


9.44 ±0.20 



Table 6. Physical parameters of f3 Vir. Estimates for density, mass and log g 
are given for two different values of Au (see text). 



Parameter 


Value 


Uncertainty (per cent) 


c'LD (mas) 




1 9 
1 .z 


J yi-u wm ) 


y.4t4 it u.zu 


Z. 1 


TTp (mas) 


91.50 ±0.22" 


0.24 


RCRq) 


1.703 ±0.022 


1.3 


TeS (K) 


6059 ± 49 


0.8 




3.51 ±0.08 


2.1 


Au (fiBz) 


72.07 ±0.10 


0.14 


pig cm-3) 


0.4028 ± 0.0081 


2.0 


M(Mq) 


1.413 ±0.061 


4.3 


log{g/cms-2) 


4.125 ±0.010 


2.4 (in g) 


Av (/iHz) 


70.5 




p(g cm-3) 


0.3851 ±0.0077 


2.0 


M(M0) 


1.351 ±0.058 


4.3 


log(5'/cm s~-^) 


4.106 ±0.010 


2.4 (in g) 



" IvanLeeuweiiltoOTl) 



4.1 Radius 



The stellar radius can be determined using our limb-darkened an- 
gular diameter and the Hipparcos parallax. 



R = 6w 



C 
27r, 



(3) 



where C is the conversion from parsecs to metres, 0ld is in radi- 
ans and TTp is in arc sec. Note that the re vised Hipparcos parallax of 
91.50 ± 0.22 mas jvan Leeuwenll200"7h has a subst antially smaller 
uncertainty than the original value of 91.74± 0.77 jPerrvman et al.l 
[1997). We obtain a stellar radius for /3 Vir of 1.703 ± 0.022 Rq. 
This is in excellent agreement with 1.706 ± 0.037Rq produced 
from surface brightne ss relations and the Hipparcos parallax by 
iThevenin et al. I( l2006l) . but is more precise and is based on a direct 
measurement. 



over the range 2.2-20 /im shows that it has negligible effect on the 
bolometric flux (<0.04 per cent) or effective temperature (<1 K). 

The resulting bolometric flux is (9.44±0.20) x 10~"' Wm"^ 
where the estimated uncertainty takes into account the fact that 
the uncertainties in the integrated fluxes for the wavelength ranges 
0.33-0.86 /im, 0.86-2.2 ^m and 2.2-20 /im are only partially in- 
dependent. Our value is consistent with two previous determina- 
tions, namely 9.408 x 10~^° Wm~^ bvlAlonso et alj lll995h and 
9.59 X 10"^" Wm"^ bv lBlackwell & Lvnas-Gravl ( Il998h" 



4 STELLAR PARAMETERS 

The observable quantities, limb-darkened angular diameter (^ld), 
bolometric flux (/) received from the star, parallax (TTp) and large 
separation {Av) combine to produce experimental constraints on 
the stellar radius {R), effective temperature (Tefj), luminosity (L), 
mean density (p), mass (Af) and surface gravity (log g). Values for 
these observable quantities are given in Table[6l along with the stel- 
lar parameters that are calculated in the remainder of this Section. 



4.2 Effective Temperature 

The combination of the stellar bolometric flux and angular diameter 
yields an empirical effective temperature: 



4/ 



1/4 



(4) 



where o is the Stefan-Boltzmann constant. Combining the new 
value for the bolometric flux (Section[3]l with the limb-darkened an- 
gular diameter (Section|2j gives the effective temperature for (3 Vir 
as 6059 ± 49 K. 

The new value for the effective temperature can be com- 
pared with the previous estimates listed in Table |3] There is 
excellent agree ment with the line-depth ratio deter mination by 
iKovtyukh etliTl ilQQj) and the determ inations by IChen et al] 
l 200lh and lAllende Prieto et al.l ( |2004|)_ from the Stromgren 
photometry using the calibratio ns by lAlonso et iD ( Il996h and 
lAUende Prieto & Lambert d 19991) . The methods are completely in- 
dependent, with only the new value presented here being based on 
a direct measurement of the angular diameter of the star. 

The values fo r t he bo lometric flux of 9.408 x 10"^" Wm~^ 
by lAlonso et al.l (Il995| ) and 9.59 x 10"^° Wm"^ by 
iBlackwell & Lvnas-GravT l l 19981) led those authors to effec- 
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live temperatures of 6088 K and 6124K respectively. However, 
if tlieir bolometric fluxes are combined with our value for the 
limb-darkened angular diameter they give effective temperatures 
of 6054 K and 6083 K respectively, both of which are consistent 
with the value for the effective temperature presented here. 



4.3 Luminosity 

The luminosity of a star can be calculated via 



(5) 



where C is the conversion from parsecs to metres and TTp is 
in arcsec. Substituting our value for / gi ves L = (1 . 349 r fc 
0.029) X 10^** W. Following the work of iBahcall et alj ( 1200 ll) . 
we have adopted Lq = (3.842 ± 0.015) x 10^'^ W and so get 
L = 3.51 ±0.08 Lq. 

A recent determinatio n of 13 Vir's luminosity by 
lEggenberger & Carried ([2006^, based on applying a bolometric 
correction to the absolute visual magnitude, gave 3.51 ± 0.10 Lq, 
which agrees with our value. 



4.4 Mean Density from Asteroseismology 

ICarrier et alj ( |2005|) collected high-precision velocity measure- 
ments of /3 Vir over eleven nights with the CORALIE spectro- 
graph. The Fourier spectrum showed a clear power excess in a 
broad frequency range centred at 1.5 mHz, confirming the pres - 
ence of solar -like oscillations first reported bv lMartic et alj j2004h . 
ICarrier et al J found that the autocorrelation of their power spectrum 
showed peaks at 70.5, 72 and 74 /iHz, and they identified the sec- 
ond of these as the most likely value for the large frequency sep- 
aration. After extrac ting oscillation frequencies and fitting to the 
asymptotic relation ( lTassoullll980h . they reported a value for the 
large separation of 72.07 ± 0.10 /^Hz. 

To a good approximation, the mean stellar density can be cal- 
culate d from the observed large frequency separation (e.g. lUkichl 
[198^ : 



(6) 



and adopting solar v alues of = 
= 134.8 /iHz tKieldsen et al.( |2008|) . 



1.408 gem""* and 
we obtain p — 



Using the large separation for f3 Vir reported by ICarrier et al.l 
l l2005h , and 

= : 

0.4028 gem"''. 

It is important to note that equation (|6j is an approxima- 
tion. In particular, Au in a given star varies systematically with 
both frequency and with the angular degree of the modes. A more 
sophisticated approach involves comparing the observed oscilla- 
tion frequencies with those calculated fr om an appropriate stel- 
lar m odel (this was done for /? Vir by lEggenberger & Carrieil 
I2OO61 ). Unfortunately, due to the difficulty in modelling convec- 
tion in the surface layers of stars, model calculations do not 
exac tly reproduce observed oscillation frequencies, even for the 
Sun jChristensen-Dalspard et alJI 19881 : iDziembowski et ai]| 19881 : 
[Rosenthal et al.lll999l: lui et alj|2002l)~ This discrepancy increases 
with frequency, which means that the large separation is incorrectly 
predicted by the mod el calculations. For example, as discussed by 
iKieldsenet ai]l l2008l) , the best models of the Sun have a large sep- 
aration that is about 1 /xHz gr eater than the observe d value. 

To address this problem. lKieldsen et aljjioo^ have proposed 



an empirical method for correcting the frequencies of stellar mod- 
els. They applied this method to the stars a Cen A, a Cen B and 
j3 Hyi and obtained very accurate estimates of the mean stellar den- 
sities (better than 0. 5 per cent). The same approach was used by 
iTeixeira et al] bOOSl) to measure the density of r Cet with similar 
precision. 

We have attempted to apply t his method to /3 Vir , using the 
oscillation frequencies published bv lCarrieret ai] l l2005h . We com- 
puted theoretical models using the Aarhus stellar evolution code 
(ASTEC, Christensen-Dalsgaard 2008a), and oscillation frequen- 
cies using the Aarhus adiab atic oscillation package (ADIPLS, 
IChristensen-Dalsgaard"2008b'). We followed the method described 
by Kjeldsen et al. (2008), which involves fitting both Av and the 
absolute frequencies of the radial modes (i.e., those having degree 
I — 0). The method makes use of the fact that the frequency offset 
between observations and models should tend to zero with decreas- 
ing frequency. However, we were not able to achieve a fit with mod- 
els whose fundamental parameters (luminosity and radius) agreed 
with the measured values. In other word s, the frequencies of the ra- 
dial modes listed bv lCarrier et al. I( l2005h do not appear to be consis- 
tent with models. Indeed, for /3 Vir Eggenberger & Carried ( |2006|) 
remarked on an offset of more than 20 fiHz between observed and 
calculated frequencies, which they noted to be substantially larger 
than the corresponding offset for the Sun. Note that our conclu- 
sion that the published frequencies are inconsistent with models 
is not sensitive to th e input physics of the evolutionary models. 
iKieldsen et al. |l|200i) demonstrated that their fitting process works 
for models published by different authors using a range of model 
codes. 

Two effects could have introduced a systematic offset to the 
observed frequencies. Firstly, an offset of ±11.6pHz could oc- 
cur because of the cycle-per-day aliases in the power spectrum, 
which are strong in single-site observations such as t hose obtained 
for (3 Vir. Indeed, many of the frequencies listed by ICarrier et al.l 
their Table 2) were shifted by 11.6/^Hz in one direction 
or the other in order to give a good fit to the p-mode spectrum. 
This includes the highest peak in the power spectrum, from which 
1 1.6 fiRz was subtracted. 

A second possibility is that the modes were misidentified, so 
that modes of odd and even degree were interchanged. That is, 
modes identified as having degree I = actually have degree I — 1, 
and vice versa. Such a misidentification is plausible and, in the con- 
text of the method of Kjeldsen et al., would be equivalent to sub- 
tracting an offset of ^Aiy ~ 36 /iHz from the frequencies of all 
modes. 

We found that neither of the above corrections was able to give 
a consistent result, but we did find a good agreement between mod- 
els and observations if we applied both together. That is, adding 
11.6/xHz and subtracting 36 fiRz from the frequencies identified 
as radial modes gave a sensible result. However, if this scenario is 
correct then the whole mode identification process would be called 
into question, and it would be risky to rely on the tabulated frequen- 
cies to establish a precise density. 

Until more oscillation data become available, we therefore fall 
back on the density that we de termined from equation Looking 
at the results for other stars dKieldsen et al.l I2OO8I : ITeixeira et al.l 
l2009l) . we estimate that this relation gives the mean density with 
an uncertainty of about 2 per cent. We have therefore adopted this 
value for the uncertainty in the density, as shown in Table |6| 
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4.5 Mass and Evolutionary State 

Combining our radius measurement witii the mean density in- 
ferred from asteroseismology, we calculate the mass to be 
1.413 ± 0.061 Mq. This is slightly higher than previous estimates 
based on fitting models in the H-R diagram: 1.36 =b 0.09 Mq 
jAllende Prieto & La mbert 199 9|), 1.32 Mq Ichen et al.l200lh and 
1.34 ± 0.10 Mq fcambert & Reddvll2004h . The difference is not 
statistically significant, but we should point out that adopting a 
value of Ai/ = 70.5 ^Hz, which is another of the three peaks in 
the autocorrelation of the observed power spectrum identified by 
ICarrier et alj ( l2005h , gives a mass of 1.35 Mq. This value, which 
is also shown in Table|6l is in better agree ment with the location of 
the sta r in the H-R diagram. The models bv lEggenberger & Carrieil 
( 12006 *) suggest that f3 Vir either has a mass of ~ 1 .3 Mq and is to- 
wards the end of its main-sequence lifetime, or else it has mass of 
~1.2Mq and is in the post main-sequence phase. Assuming that 
one of 70.5 or 72. 1 fiKz is the correct large spacing, then the post 
main-sequence model appears to be excluded, indicating that /3 Vir 
is still on the main sequence. Clearly, establishing the age of the 
star is also affected by the uncertainy in the mass. A more precise 
determination of its mass and evolutionary state will require further 
observations of its oscillations, preferably from multiple sites. 

Substituting the stellar mass with the mean density and vol- 
ume, the standard surface gravity relation becomes 

g = ^GttpR, (V) 

where G is the universal constant of gravitation. In Table[6]we show 
log g for the two values of Av. Both determinations are consistent 
with the published values (see Table [3](, which is not surprising 
given that spectroscopic determinations of log g are not very pre- 
cise. 



5 CONCLUSION 

We have presented the first angular diameter measurement of the 
F9 V star /? Vir. In combination with the revised values for the bolo- 
metric flux and parallax, this angular diameter has experimentally 
constrained the stellar radius and effective temperature. The radius 
measurement, combined with the mean density determined from 
asteroseismology, allows an estimate of /? Vir's mass and surface 
gravity. 

The constraints on R, Ted, M, log g and L that we present will 
be invaluable in the future to further test theoretical models of B Vir 
and its oscillations. As stressed bv ISrown & Gillilan j jl994h . for 
example, oscillation frequencies are of most importance for testing 
evolution theories when the other fundamental stellar properties are 
well-constrained. 
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